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Abstract 

The gradient scheme framework is based on a small number of properties and encompasses a large number of numerical 
methods for diffusion models. We recall these properties and develop some new generic tools associated with the gradient 
scheme framework. These tools enable us to prove that classical schemes are indeed gradient schemes, and allow us to 
perform a complete and generic study of the well-known (but rarely well-studied) mass lumping process. They also allow 
an easy check of the mathematical properties of new schemes, by developing a generic process for eliminating unknowns 
via barycentric condensation, and by designing a concept of discrete functional analysis toolbox for schemes based on 
polytopal meshes. 
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1 Introduction 

A wide variety of schemes have been developed in the last few years for the numerical simulation of anisotropic diffusion 
equations on general meshes, see [23, 47, 42] and references therein. The rigorous analysis of these methods is crucial to 
ensure their robustness and convergence, and to avoid the pitfalls of methods seemingly well-defined but not converging 
to the proper model [43, Chapter III, §3.2]. The necessity to conduct this analysis for each method and each model has 
given rise to a number of general ideas which are re-used from one study to the other; a set of rather informal techniques 
has thus emerged over the years. 

It is tempting to push further this idea of “set of informal similar techniques”, to try and make it a formal mathematical 
theory. This boils down to finding common factors in the studies for all pairs {method,model), and to extract the core 
properties that ensure the stability and convergence of numerical methods for a variety of models. Identifying these core 
properties greatly reduces the work, which then amounts to two tasks: 

Task (1): establish that a given numerical method satisfies the said properties; 

Task (2): prove that these properties ensure the convergence of a method for all considered models. 

Thus, the number of convergence studies is reduced from [Ca,Td{methods)x Card(modeZs)] - which corresponds to one per 
pair {method,model) - to [Card(met/iods) -I- Card(mode/s)]. Ca,rd{methods) studies are needed to prove that each method 
satisfies the core properties, and Card( mode/s) studies are required to prove that an abstract method that satisfies the 
core properties is convergent for each model. 

Attempts at designing rigorous theories of unified convergence analysis for families of numerical methods are not new, see 
e.g. [16, 21, 9, 11] for finite element, discontinuous Galerkin methods and compatible discretisation operators. Recently, 
the gradient scheme framework was developed [36, 29]. Not only does this framework provide a unifying framework for a 
number of methods (Task 1) - conforming and non-conforming finite elements, finite volumes, mimetic finite differences, 
... - but it also enables complete convergence analyses for a wide variety of models of 2nd order diffusion PDEs (Task 2) 
- linear, non-linear, non-local, degenerate, etc. [29, 39, 31, 26, 25, 2, 27, 13, 40, 15] - through the verification of a very 
small number of properties (3 for linear models, 4 or 5 for non-linear models). 

The purpose of this article is to bring gradient schemes one step further towards a unification theory. Indeed, we develop 
a set of generic tools that make Task (1) extremely simple for a great variety of methods. In other words, using these 
tools we can produce short but complete proofs that several numerical methods for 2nd order diffusion problems are 
gradient schemes. 
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The paper is organised as follows. In Section 2, we present the gradient scheme framework. This framework is based on 
the notion of gradient discretisation, which defines discrete spaces and operators, and on five core properties, presented in 
Subsection 2.1: coercivity, consistency, limit-conformity, compactness, and piecewise constant reconstruction. A gradient 
scheme is a gradient discretisation applied to a given diffnsion model, consisting in a set of second order partial differential 
equations and bonndary conditions. Depending on the considered model, a gradient discretisation must satisfy three, 
four or five of these core properties to give rise to a convergent gradient scheme. In sections following 2.1, we develop 
generic notions that are useful to establish that particular methods fit into the framework. More precisely, in Subsection 
2.2 we introduce the concept of local linearly exact gradients, and we show that it implies one of the core properties - the 
consistency of gradient discretisations. Subsection 2.3 deals with the barycentric condensation of gradient discretisations, 
which is a classical way to eliminate degrees of freedom. The gradient scheme framework enables us, in Subsection 2.4, 
to rigorously define the well-known technique of mass lumping, and to show that this technique does not affect the 
convergence of a given scheme. In Subsection 2.5 we provide an analysis toolbox for schemes based on polytopal meshes, 
and we introduce the novel notion of control of a gradient discretisation by this toolbox. This notion enables us to 
establish three of the main properties (coercivity, limit-conformity and compactness) and therefore completes the notion 
of local linearly exact gradient discretisations. 

In Section 3, we show that all methods in the following list are gradient schemes and satisfy four of the five core 
properties (coercivity, consistency, limit-conformity, compactness): conforming and non conforming finite elements, RTfe 
mixed finite elements, multi-point flux approximation MPFA-O schemes, discrete duality finite volume (DDFV) schemes, 
hybrid mimetic mixed methods (HMM), nodal mimetic finite difference (nMFD) methods, vertex average gradient (VAG) 
methods. For these methods, the fifth property (piecewise constant reconstruction) is either satisfied by definition, or 
can be satisfied by a mass-lumped version in the sense of Subsection 2.4. The mass-lumped versions are only detailed 
in the important cases of the conforming and non-conforming Pi finite elements. We show that the notions of local 
linearly exact gradient discretisations, and of control by polytopal toolboxes, apply to most of the considered methods, 
and therefore provide very quick proofs that these methods satisfy the consistency, coercivity, limit-conformity and 
compactness properties. Some of the schemes have already been more or less formally shown to be gradient schemes in 
[36, 29], but the proofs provided here thanks to the new generic tools developed in Section 2 are much more efficient and 
elegant than in previous works, and can be easily extended to other schemes. 

A short conclusion is provided in Section 4. 

2 Gradient discretisations: definitions and analysis tools 

For simplicity we restrict ourselves to homogeneous Dirichlet boundary conditions; all other classical boundary conditions 
(non-homogeneous Dirichlet, Neumann, Fourier or mixed) can be dealt with seamlessly in the gradient schemes framework 
[27]. The principle of gradient schemes is to write the weak formulation of the PDF by replacing all continuous spaces 
and operators by discrete analogs. These discrete objects are described in a gradient discretisation. Once a gradient 
discretisation is defined, its application to a given problem then leads to a gradient scheme. 

For linear models, the convergence of gradient schemes is obtained via error estimates based on the consistency and 
limit-conformity measures St> and W-d- For non-linear models, whose solutions may lack regularity or even be non¬ 
unique, error estimates may not always be obtained; however, convergence of approximate solutions can be obtained 
via compactness techniques such as those developed in the finite volume framework [32, 35, 23]. Even though they do 
not yield an explicit rate of convergence, these compactness techniques provide strong convergence results - such as 
uniform-in-time convergence [25] - under assumptions that are compatible with field applications (discontinuous data, 
fully non-linear models, etc.). 

2.1 Definitions 

Definition 2.1 (Gradient discretisation for homogeneous Dirichlet boundary conditions): Letp G (l,oo) and let 17 be a bounded 
open subset of R'*, where d £ N\ {0} is the space dimension. The triplet V — (A-d.o, n^j, V-d) is a gradient discretisation 
for problems posed on 17 with homogeneous Dirichlet boundary conditions if: 

1. Xt),o is a finite dimensional space encoding the degrees of freedom (and accounting for the homogeneous Dirichlet 
boundary conditions), 

2. Ht) : Xt>,o —t L^iyi) is a linear mapping reconstructing a function in U’iyi) from the degrees of freedom, 

3. V-D : X-Dfi —t U'{Q.Y is a linear mapping defining a discrete gradient from the degrees of freedom, 

4- IjV-D • ]|LP(n)ti is a norm on Axi,o. 

Here are the three properties a gradient discretisation needs to satisfy to enable the analysis of the corresponding gradient 
scheme on linear problems: 

• The coercivity ensures uniform discrete Poincare inequalities for the family of gradient discretisations; this is 
essential to obtain a priori estimates on the solutions to gradient schemes. 
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• The consistency states that the family of gradient discretisations “covers” the whole energy space of the model 
(e.g. Hq{Q.) for the linear equation (2.1)). 

• The limit-conformity ensures that the family of gradient and function reconstructions asymptotically satisfies the 
Stokes formula. 


Definition 2.2 (Coercivity): Let be a gradient discretisation in the sense of Definition 2.1 and let C-d be the norm of 
the linear mapping IIt> defined by 


Ct> 


lini5M||iP(n) 

max -77=-n- 


A sequence (T>m)meN of gradient discretisations in the sense of Definition 2.1 is said to be coercive if there exists Cp £ R+ 
such that C-Vm ^ Ctp for all m £ N. 


Definition 2.3 (Consistency): Let be a gradient discretisation in the sense of Definition 2.1 and let St> '■ ITg’^(D) —> 
[0, +oo) be defined by 

^(p € Wo’^{Q), ST>i(p) = min (HII-du - (p||pp(n) + || Vpw - Vv9||pp(n)d) . 

uGXd^o 

A sequence (Dm)meN of gradient discretisations in the sense of Definition 2.1 is said to be consistent if for all ip £ ITg’^(D) 
we have limm^oo >S'p„(v5) = 0. 

Definition 2.4 (Limit-conformity): Let D be a gradient discretisation in the sense of Definition 2.1. We set p' = the 
dual exponent of p, and (D) = {ip £ U’ (D)‘*, div(y5 £ U' (D)}, and we define 

(D), Wt}{(p) = sup - f iy-Duix) ■ ip{x) +Il'Du{x)Ahfifi{x))Ax . 

“eXp,_o\{o} II V■DM||pp(n)ci Jet 

A sequence (Dm)meN of gradient discretisations is said to be limit-conforming if for all tp £ (O) we have 

lim^^oo W-D„((^) = 0. 


To give an idea of how a gradient discretisation gives a converging gradient scheme for diffusion equations, let us consider 
the case of a linear elliptic equation 

r -div(AVM) = / in D, , . 

\ u = 0 on do., ^ ’ 

where A : D i—>■ A1d(R) is a measurable bounded and uniformly elliptic matrix-valued function such that A{x) is 
symmetric for a.e. x £ Q,, and / £ L^{Q). The solution to problem (2.1) is understood in the weak sense: 

Find u £ Hq{Q) such that, for all v £ Ro(D), / A(x)\/u{x) ■ Vu(a;)da: = / f{x)v{x)dx. (2-2) 

Jn Jn 

If D is a gradient discretisation with p = 2, then the corresponding gradient scheme for (2.1) consists in writing 


Find u £ Xt>,o such that, for all v £ X-pp, 


A{x)Xpu{x) ■ Xpv{x)dx = / f{x)Tlpv{x)dx. 


I 


(2.3) 


As seen here, (2.3) consists in replacing in (2.1) the continuous space Hq{Q) and the continuous gradient and function 
by their discrete reconstruction from D. Reference [36] proves the following error estimate between the solution to (2.2) 
and its gradient scheme approximation (2.3): 


||Vm — Vi)M||^2(-f2)‘i + l|w — n-pull < Cl [Wp{AXu) Sp{u)] , (2.4) 

where Ci only depends on A and an upper bound of Cp in Definition 2.2. This shows that if a sequence (Dm)meN of 
gradient discretisations is coercive, consistent and limit-conforming, and if (um)meN is a sequence of solutions to the 
corresponding gradient schemes for (2.1), then Ilp^Um —>■ u in L^(D) and Xp^Um —>■ Xu in The study of a 

scheme for (2.1) then amounts to finding a gradient discretisation T> such that the scheme can be written under the form 
(2.3), and to proving that sequences of such gradient discretisations satisfy the above described properties. Establishing 
the consistency and limit-conformity usually consists in obtaining estimates on Sp and Wp that give explicit rates of 
convergence in (2.4). 

Dealing with non-linear problems might additionally require one or both of the following properties. 

• The compactness is used to deal with low-order non-linearities - e.g. in semi- or quasi-linear equations. 

• The piecewise constant reconstruction corresponds to mass-lumping and is required to manage certain monotone 
non-linearities, or non-linearities on the time derivative as in Richards’ model. 
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Definition 2.5 (Compactness): A sequence (Vm)meN of gradient discretisations in the sense of Definition 2.1 is said to be 
compact if, for any sequence (um)meN such that Um € for m £ N and (||V'D„Mm||iP(n)<i)m6N is bounded, the 

sequence (nj>^Um)meN is relatively compact in L’’{Q). 

Definition 2.6 (Piecewise constant reconstruction): Let D = (X-d^o, II-d, V-d) be a gradient discretisation in the sense of 
Definition 2.1. The linear mapping II-d : Xt>,o —t L^(D) is a piecewise constant reconstruction if there exists a 
finite set B, a basis {ei)i£B of Xr>_o and a family of (possibly empty) disjoint subsets (Vi)igB of D such that for all 
u — Uid € Xt>,o we have Tl-vu = tiiXVi, where XVi is the characteristic function of Vi. 

Remark 2.7: Piecewise constant reconstructions generally use as set B the set I of geometrical entities attached to the 
degrees of freedom of the method (see Definition 2.11); in this case (ei)igs is the canonical basis of X-pfl- Note that it is 
possible, starting from a generic gradient discretisation B, to replace the original reconstruction \It> by a reconstruction 
that is piecewise constant; the new gradient discretisation thus obtained is called a mass-lumped version of T> (see Section 
2.4). 


As an illustration of the use of the importance of these properties for nonlinear problems, let us consider the following 
semi-linear modification of (2.1): 

-div(AVM) -(- I3{u) = / in D, 

M = 0 on dVL, ' ' 

for some function /3 such that /3(s)s > 0 for all s £ R. The gradient discretisation of this problem is pretty straightforward: 
find u £ X-pfi such that, for all v £ Xt>,o, 


A(x)Xdu{x) ■ X'dv{x)Ax -\- / p{Il-r>u(x))Il-r>v{x)dx = / f{x)Il-r>v{x)dx. 

! Jn Jn 


( 2 . 6 ) 


The compactness property implies that an estimate on a sequence of discrete gradient (V-Pm'>^rn)meN will yield relative 
compactness of the corresponding sequence of reconstructions (n-D^Um)meN, thus enabling a passing to the limit in the 
nonlinearity /3. 

The formulation (2.6) ensures a priori estimates on the solution since, taking v = u, the term j3(JlT>u)YlT)U is non¬ 
negative. However, in practical implementation, computing /3(nr>u(x))nT>v(x)dx might be problematic; even if Yl-pu 
and YIt>v are polynomials on each cell of a mesh (as in finite element schemes), j3{JlT>u)\lT>v is locally not polynomial 
and no exact quadrature rule might exist to compute its integral. An alternative scheme consists in replacing, in (2.6), 
the term 

f l3(JlT>u{x))IlT)v{x)dx (2.7) 


with 


YlT>l3{u){x)Ii'Dv{x)dx, 


( 2 . 8 ) 


where /?(«) G Xt>,o is defined degree-of-freedom per degree-of-freedom, that is {P{u))i = p{ui) for all i £ B with the 
notations of Definition 2.6. For finite element methods, (2.8) consists in integrating polynomials on cells, and exact 
quadrature rules can be used. However, this alternative scheme does not ensure a priori estimates on the solution since, 
with the choice v = u, the term might be negative in some part of the domain. Hence, we have to choose 

between unconditional stability (o priori estimates) with (2.7), or a scheme that is practical to implement with (2.8). 
One of the interests of piecewise constant reconstructions is to solve this apparent contradiction. If Ut> is a piecewise 
constant reconstruction then 

ni,(/3(u)) = I3{I1t,u). (2.9) 

Hence (2.7) and (2.8) are identical, and both stability and computational practicality are satisfied. The second interest 
of piecewise constant resconstructions can be found in the analysis of time-dependent problems. The discretisation of 
dtu leads to a term of the form 

f nBu"+\x)-nT,u"{x)^ 


-Ii'Dv(x)dx. 


)i£B 


( 2 . 10 ) 

of in (2.10) is 


Jn ^ 

If H-d is a piecewise constant reconstruction, the mass matrix multiplying the coordinates (m"^^) 
diagonal, and its inversion is therefore trivial. 

As shown in [25, 31, 39], piecewise constant reconstructions ensure the stability and convergence of gradient schemes for 
a variety of non-linear elliptic and parabolic equations. 


Remark 2.8: 1. The consistency, limit-conformity and compactness of gradient discretisations may be defined in other 

equivalent ways [27]. Moreover, the consistency of a sequence of gradient discretisations only needs to be checked 
for ^ in a dense set of the domain of St> (e.g. C“(D)). The limit-conformity of a coercive sequence of gradient 
discretisations only needs to be checked for in a dense set of the domain of W-v (e.g. C“ (R‘*)‘*, which is indeed 
dense in (D) when D is locally star-shaped, which is the case if Q is polytopal). Finally, the compactness 

of a sequence of gradient discretisations implies its coercivity. 

2. Gradient discretisations for time-dependent problems can be easily deduced from the gradient discretisations for 
steady-state problems [29, 27]. 
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2.2 Local linearly exact gradients 

Most numerical methods for diffusion eqnations are based, explicitly or implicitly, on local linearly exact reconstructed 
gradients. The following definition gives a precise meaning to this. 


Definition 2.9 (Linearly exact gradient reconstructions): Let U he a bonnded set of R'*, let / be a finite set and let S = 
{xi)i^i be a family of points of R'*. A linear mapping Q : R^ i—>■ L°°{U)^ is a linearly exact gradient reconstruction npon 
S if, for any affine fnnction L : R'* —>■ R, if ^ = {L{xi))i^i then = VL on U. The norm of Q is defined by 


lien 


diam([/) max 
5eRr\{0} 


maxig/ |^i| 


( 2 . 11 ) 


As expected, linearly exact gradient reconstrnctions enjoy nice approximation properties when compnted from inter- 
polants of smooth fnnctions. 

Lemma 2.10 (Estimate for linearly exact gradient reconstructions): Let [/ be a bounded set of R'*, let S — {xi)i^i C R'*, 
and let Q : R^ i— L°°{U)^ be a linearly exact gradient reconstrnction upon S in the sense of Definition 2.9. Let 
ip G 1T^’°“(R‘*) and define v G R^ by Vi = <p{xi) for any i £ I. Then 

\gv - Vy| ^ (^ + ^II^IU j diam([/)||v3||w2.=o(R<i) a.e. on t/. 


Proof Take xu G U and let L{x) = ip(xu) + \/p{xu) ■ {x — xu) be the first order Taylor expansion of ip aronnd xu- 
Let ^ = {L(xi))i^i. By linear exactness of Q we have = VL = \/p{xu) on U. Hence, 

I^C - Vv9| < diam(t/)||(p||w 2 ,cx,(R,d) onU. (2.12) 

For any f G / we have (v — ^)i = p{xi) — L{xi) = p{xi) — p{xu) — '^p{xu) ■ {xi — xu). Since \xi — xu\ < dist(a:i, U) + 
diam(t/), we get |(v — ^)i| < |(dist(a:i, U) + diam(f7))^||g)||j.(^2,oo(H£i). The linearity of Q and the definition of its norm 
therefore imply, for a.e. x £ U, 

\gv{x) - g5(a:)| = - e(a;)| < \ C) + diam(C/)i 

< l|l5|UdiR.„(C) + l)’ 

Combined with (2.12), this completes the proof of the lemma. □ 

The consistency of gradient discretisations based on linearly exact gradient reconstrnctions follows. Let us first give the 
the definition of snch gradient discretisations. 


Definition 2.11 (LLE gradient discretisation): The triplet V = (Xc^o, H-d, Vr>) is an LLE (for “local linearly exact”) gra¬ 
dient discretisation if there exists a finite partition W of fi, a set / of geometrical entities attached to the degrees of 
freedom (dof), a finite family of approximation points S — (xi)i^i C R'* and, for any U £14, a. subset lu <£ I such that: 

1. Xt>,o ~ R^^ X where the set / is partitioned into In (interior geometrical entities attached to the dof) and 

Ian (bonndary geometrical entities attached to the dof). 

2. There exists a family {ai)i^i such that, for all i £ I, ai £ L°°(fl) and 

(a) Vf G I, yU £l4, li i ^ lu then Oi = 0 on U, 

(b) for a.e. x £ fl, ai(x) = 1 and Vv £ Xx>,o, hi'Dv{x) = ai{x)vi. (2.13) 

i£l i£l 

3. There exists a family {gu)ueu snch that, for all U £ 14, Qjj : R^'^ i—>■ is a linearly exact gradient 

reconstrnction upon {xi)i(zjj^, in the sense of Definition 2.9, and V-dv = gu{ivi)ieiu) on U, for all v £ X-v^. 

4. IjVr, • II LP(n)a is a norm on Xt},o- 

In that case, we define the LLE regularity of T> by 

regLLEp) = max f WGuWo^ + max ) -L esssup^ |ai(a;)|. (2.14) 

ueu \ leiu aiamfty j J ^sen 

Remark 2.12: As implied by the terminology, an LLE gradient discretisation is a gradient discretisation in the sense of 
Definition 2.1. Note that the existence of i,j £ I with i ^ j and Xi = Xj is not excluded (see, e.g.. Section 3.4). 
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Remark 2.13: We do not request HdV to be linearly exact (oi is not necessarily affine in each U)', this reconstruction just 
needs to be computable from local degrees of freedom, and exact on interpolants of constant functions. This enables us 
to consider mass-lumped gradient discretisations. 

In a number of cases, estimating s--®- * G is trivial. For example, if for a.e. x £0., there is exactly one 

i £ I such that ai{x) — 1 and Qj(x) = 0 for all j £ we get ~ 1 s--®- (then ID has a piecewise constant 

reconstruction and the set B defined in Definition 2.6 is identical to I). Another example is the case where, for a.e. 
X £ U, B-uvix) is a convex combination of the dof {vi)i£i^ (which is the case, e.g., if Tljyv is linear on U, Vi = IlT>v{xi) 
and (xi)i^ijj are extremal points of U)\ then Oi > 0 for all f G / and l'^*(®)l = 1- 

Proposition 2.14 (LLE gradient discretisations are consistent): Let (D>m)mgN be a sequence of LLE gradient discretisations, 
associated for any m G N to a partition Um- If (regLLg(lDm))mgN is bounded and if max diam(t/) —>■ 0 as m —>■ oo, then 

(D>m)mgN is consistent in the sense of Definition 2.3. 

Proof Let tp £ C“(fl) and let w"* = (</9(a;’"))igj™ G where Sm = (a;™)igj™ is the family of approximation 

points of T>m- Owing to Lemma 2.10 we have, for U £ Um and a.e. x £ U, 

- Vp{x)\ = - Vp{x)\ 

- + ^''®SLLE(®m)(regLLE(®™) + 1)^^ diam(t/)||(p||,v2,cx=(R,d). (2.15) 

Let us now evaluate — p\. Since any (a;7*)igj™ is within distance regLLE(D>m)diam(f7) of U, for all i £ I™ and 

all X £U we have \v™ — p{x)\ < (1 -I- r®gLLE(®™))diam(t/)||(p||,.(..i,oo(pi£i). By (2.13), we infer that for a.e. x £ U 


|II-D„i’™(a:) - ‘Pix)\ 


aT{x){vr - p{x)) 

iei^ 


< \(b{x)\ 


I m 

sup \Vi 


lfi{x)\ 


< regLLE(®m)(l + i'®gLLE(®m))diam(t/)||(p||,vi,oo(R,d). 


(2.16) 


Estimates (2.15) and (2.16) and the assumptions on (IDm)meN show that —>■ Xip in L°°(fl)‘* and —>■ (f 

in L°°{Q) as m —>■ 00 . Remark 2.8 then concludes the proof. □ 


2.3 Barycentric elimination of degrees of freedom 

The construction of a given scheme often requires several interpolation points. However, some of these points can be 
eliminated afterwards to reduce the computational cost. A classical way to perform this reduction of degrees of freedom 
is through barycentric combinations, by replacing certain unknowns with averages of other unknowns. We describe here 
a way to perform this reduction in the general context of LLE gradient discretisations, while preserving the required 
properties (coercivity, consistency, limit-conformity and compactness). 


Definition 2.15 (Barycentric condensation of an LLE gradient discretisation): Let D be an LLE gradient discretisation. We 
denote by S' = {xi)i^i C R'* the family of approximation points of V and by U its partition. A gradient discretisation 
is a barycentric condensation of B if there exists C I and, for all i £ a set Hi C and real numbers 

satisfying 

= 1 and (2-17) 

jeHi jeH, 

such that 

• Ian C 

• XpB.^o = X 

• For all V £ XpBA Q we have lipBAV = II-dE and Vj^baU = VpV, where V £ Xpfi = R^“ x {0}^®^ is defined by 


Vi G /, Vi ^ 


Vi 

E, 


if i G 

(We note that V is indeed in Xpfi since Ian C and Ui = 0 if i G Ian.) 
We define the regularity of the barycentric condensation B^'' by 


(2.18) 


r®gBA(2^‘'*) = 1 + max 






distfajy, Xi 

max max-^ 

UGU \ i€iu jGHi diam((7) 
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It is clear that the above defined barycentric condensation O®* is a gradient discretisation. Indeed, if Vj^baV = 0 on 
then V-dV = 0 on and thus Vi = 0 for all f £ S' (since V is a gradient discretisation and therefore HV-d • ||LP(n)<i is a 
norm on This shows that Vi = 0 for all i £ and thus that HVj^ba • ||_LP(n)d is a norm on Xj^ba q. 

Remark 2.16 (Localness of the barycentric elimination): Bounding the last term in reg0,^(I2“*) consists in requiring that, if 
i £ 1 \ is involved in the definition of Qu, then any degree of freedom j £ Hi used to eliminate the degree of freedom 
i lies within distance (!l(diam([/)) of U. This ensures that, after barycentric elimination, Qu is still computed using only 
degrees of freedom in a neighborhood of U. 

Barycentric elimination expresses some degrees of freedom by combinations that are linearly exact. As a consequence, the 
LLE property is preserved in the process, and the consistency of barycentric condensations of LLE gradient discretisations 
is ensured by Proposition 2.14. 

Lemma 2.17 (Barycentric elimination preserves the LLE property): Let V be an LLE gradient discretisation in the sense of 
Dehnition 2.9, and let be a barycentric condensation of T>. Then is an LLE gradient discretisation on the same 
partition as V, and < regB,^(I>®*) regLLE(®) + regBA(I>°*) + regLLE(®)- 

Proof Obviously, = (7®* 0 7n) U 7gn forms the geometrical entities attached to the dof of P®* since X^^ba q = 
Let 77 be the partition corresponding to T>, and let U G U. Take v £ X^^ba q and let V £ 
be defined by (2.18). We notice that, for any U G U, the values {Vi)i£i^ are computed in terms of with 

7&^ = (7r;n7-)uUe,^\,BA77.. 

We have, for x G U, 

ni,BAu(a;) = 11151^(0;) = ^ ai{x)Vi = ^ ai{x)vi + ^ ai(x) E = E ai{x)vi 

with 

ai(x) = ai{x) + ^ piUkix) if i £ 7^ n 7®*, 

keiu\i^^\ieHk 

ai{x) = ^ Piakix) itiGllP\Iu- 

Thanks to (2.17) and (2.13) we have 


E ai(®) = E E E PiOikix) 

i6%ri/BA \ieHi. 

= ^ ai{x) + ^ ak{x) Pi ^ E + E = 1- (2-19) 


Hence, UpBAV has the required form. The gradient {'VpBAv)\u = Qu{{Vi)i<=ijj) only depends on (ui)iG/))* thus be 

written Qu{{'>^i)i^iBA). By (2.17) the reconstruction v i—>■ E is linearly exact, that is if v interpolates the values of an affine 
mapping L at the points (xiJ^^jBA then V interpolates the same mapping L at the points {xi)i^i^. Hence, the linear 

exactness of Qu gives the linear exactness of Gu- This completes the proof that P®* is an LLE gradient discretisation. 

Let us now establish the upper bound on regLLE(®"^*)- For all i G we have \Vi\ < l/^tl l^ll — ''OSba(F’"**) rnax^gjBA 

This also holds for i G lu 7®* since regBA(T’^*) > L Hence, a.e. on U, 


9u {{vi)i^jBA^ = \gu {(yi)ieiu)\ < 


Hgi/||ooregB^(P‘^^) 

diam(7/) 


max luil 
iei^-^ 


and thus 

WGuU < liec/||ooregBA(7)“"). (2.20) 

Reproducing the reasoning in the first two equalities in (2.19) with absolute values and inequalities, we see that 


E - E i«i(®)i + E i«*(®)i E E (2.21) 

ieijp' ieiucii^^ keiu\i^‘' i&iu 

Finally, for j G Ifp we estimate by studying two cases. If j £ lu then dist(a:j,7/) < regLLB(T’)diam(7/). 

If j ^ lu then there exists i G Iu\I^^ such that j G Hi, and thus dist(a:j, cci) < regBA(®°*)diam(7/); this gives 
dist(a;j,7/) < (regBA(T’*^^) + regLLB(7?))diam([/). Combined with (2.20) and (2.21), these estimates on dist(a:j,[/) prove 
the bound on regLLE(F’°^) stated in the lemma. □ 

Barycentric condensations of LLE gradient discretisations satisfy the same properties (coercivity, consistency, com¬ 
pactness, limit-conformity) as the original gradient discretisation. The coercivity, limit-conformity and compactness 
properties result from the fact that Xj^ba ^ is (roughly) a subspace of X-d,o, and the consistency is a consequence of 
Lemma 2.17 and Proposition 2.14. 



2 Gradient discretisations: definitions and analysis tools 


Theorem 2.18 (Properties of barycentric condensations of gradient discretisations): Let (Vm)meN be a sequence of LLE gra¬ 
dient discretisations that is coercive, consistent, limit-conforming and compact in the sense of the definitions in Section 
2.1. Let Ltm be the finite partition of corresponding to ’Dm- We assume that maxusUm diam(t/) —>■ 0 as m —>■ cx 3 , 
and that (regLLg(Ilm))m 6 N is bounded. For any m G N we take a barycentric condensation of Om such that 

(regB,^(T>;)^))m 6 N is bounded. 

Then is coercive, consistent, limit-conforming and compact. 

Proof For any v £ Xjjba q, with V defined by (2.18) we have 

l|n-D^Au||LP(n) = ||ni5„^P||iP(n) < C-D„||Vi:)„P||^p(n)d = 11Vi,baIi| |^p(f2)d, 

which shows that Cjjba < C-Vm S'lid thus that (T>))^)m 6 N is coercive. To prove the compactness, we take (V 2 ,BAWm)m 6 N = 
(V'D„Kri)m 6 N bounded in L^(fl)‘*, and we use the compactness of ('Dm)meN to see that (n-D„Kp)mgN = (n 2 ,BAWm)meN 
is relatively compact in The limit conformity follows by writing 




J ■ (fi{x) + Ylj,^v{x)diYLp{x)^ dx 


l|Vi,„P|| LP(n)'>- 


(X'Drr,V(x) ■ ^pix) n-D„P(a:)div(^(:c)) drc 


which shows that Wx>ba{p) < W-Dm{p)- Finally, by Lemma 2.17 each T>^ is an LLE gradient discretisation and the 
boundedness of (regLLj;(Iim))meN and (regB,^(T>m))meN show that (regLLE(Iim))meN is bounded. Proposition 2.14 then 
gives the consistency of (r>))^)m 6 N- □ 


2.4 Mass lumping and comparison of reconstruction operators 

“Mass-lumping” is the generic name of the process applied to modify schemes that do not have a built-in piecewise 
constant reconstruction, say for instance the Pi finite element scheme (see Section 3.1). This is often done on a case-by- 
case basis, with ad hoc studies. The gradient scheme framework provides an efficient generic setting for performing this 
mass-lumping. The idea is to modify the reconstruction operator so that it becomes a piecewise constant reconstruction; 
under an assumption that is easy to verify in practice, this “mass-lumped” gradient discretisation can be compared with 
the original gradient discretisation, which ensures that all properties required for the convergence of the mass-lumped 
scheme are satisfied. 

Definition 2.19 (Mass-lumped gradient discretisation): Let D = (Xu^o, fl-D, Vn) be a gradient discretisation in the sense 
of Definition 2.1. A mass-lumped version of H is a gradient discretisation = (Xp^o, If™, Vr>) such that If™ is a 
piecewise constant reconstruction in the sense of Definition 2.6. 

In all the cases of mass-lumping considered in this paper, we show that the following theorem applies to — 2?“. 
This theorem states that, if two sequences of gradient discretisations share the same space and reconstructed gradients, 
one inherits the properties from the other provided that their reconstruction operators are close to each other (condition 
(2.22)). Moreover, it also establishes that the sufficient condition (2.22) is also necessary for the mass-lumped schemes to 
satisfy the compactness and limit-conformity properties, since these properties are satisfied by all the considered initial 
schemes. 


Theorem 2.20 (Comparison of reconstruction operators): Let (T>m)meN be a sequence of gradient discretisations in the 
sense of Definition 2.1. ForanymG N, let be a gradient discretisation defined from by O™ = (Xi 5 ^_o, , V-d„ ), 

where LIJ,^ is a linear operator from to U’iyi). 

1. We assume that there exists a sequence (cJ„i)meN such that 


lim™_>oo iAJm = 0 , and 

Vm G N , Vu G Xt)^,o , HLIo^-u — nxi„i>||LP{n) < cu m l|Vc mf |Lp{n)£i- 

If (T’m)m 6 N is coercive (resp. consistent, limit-conforming, or compact - in the sense of the definitions in Section 
2 . 1 ), then (T>m)m 6 N is also coercive (resp. consistent, limit-conforming, or compact). 

2. Reciprocally, if (T’m)meN and are both compact and limit-conforming in the sense of the definitions in 

Section 2.1, then there exists (cUm)m 6 N such that (2.22) holds. 


Proof Let us prove the first item of the theorem. We let M = sup,^gi,,f cUm, and we use the triangular inequality to 
write, from (2.22), for any v G X-d„,o, 

||niJ„'w||_LP(n) < llllp^u — \dT}^v\\LP(n) + ||IIi:)„u||iP(n) < M|| Vo„u||£P(n)d -|- ||ni5„u||iP(o). 


(2.23) 
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COERCIVITY: let US assume that (Vm)meN is coercive with constant Cp. Using (2.23) we find ||n|,^t;||/^p(n) < (M + 
C'p)||V-D„u||iP(n)d and the coercivity of follows. 

Consistency: let us assume that (iDm)m6N is consistent. Using the triangular inequality and (2.22), we write, for any 
V G and (p e 1U((’'’(!^), 

— ^\\lp{q.) + l|Vi5„u — Vv?||iP(n)rf 

< Wm||Vi5„w||ip(n)d + llUxi^u - (p|Up(n) + ||Vi5„u - V¥?||LP(n)<i 

< Wm||V(p||2,P(n)£i + W'^'Drn'^ “ V5||iP(n) + (1 + ‘^m)||V-D„tI — V(p||^p(Q)d 

< Wm||V(p||2.p(n)£i + (1 + .iU)(||ni>„jU — ipWhPin) + ||V-d„u — V(p||2.p(n)£i). 

Hence S-oi^i^p) ^ + (1 + M)ST>m{<fi) and the consistency of (V^jmeN follows from the consistency of 

(®m)meN and from lim™_>oo Wm = 0. 

Limit-conformity: let us now assume that (Om)meN is limit-conforming. By the triangular inequality and (2.22), for 
any (p € (f2), 


J {yTimVix) ■ ip{x) + nj,^u(a;)div:^(a;)jda: 


< ll<iiv<*3||^p-(s.;)W™||Vo„u||ip(n)d -h 


/ (V'Dpit’(a:) • ‘fiix) -\- ]l-D^v{x)Aw(p{x)) Ax 
Jfi 


Using (2.22), we infer that W-d^ {ip) < cJm||div(^||p;(-fj)d A-W-Pm (‘p) —t 0 as m —>• oo, and the limit conformity of (V^)meN 
is established. 

Compactness: we now assume that (Vrn)meN is compact. If (VT 7 ^Vm)meN is bounded in U’iyiY, then the compactness 
of (U’m)meN ensures that (n-D„Um)m6N is relatively compact in L^(f7). Since \\Tl'^^Vm — n-D^Um||np(n) —>■ 0 as m —>■ cxd 
by (2.22), we deduce that (n^^ii„i)m6N is relatively compact in L^(n). 

Let us now turn to the proof, by way of contradiction, of the second item. We therefore assume that (l?m)meN and 
mgN are both compact and limit-conforming, and that 


cJm := max 

peXp, ,o\{o} 


ll^'D„'i’||iP(n)d 


0 as m —>■ oo. 


(2.24) 


Then we can find eo > 0, a subsequence of (T’m,Om)meN (not denoted differently) and for each m G N an element 
Vm G such that - nxi^UmlUp(n) > eo||Vi>„Um||iP(n)<i- Since u™ yf 0, we can consider = 


I Ppi 1 1 


LP(U )<1 


which satisfies 11V-D^Wml|iP(n)d = 1 and 


IHp^Wm — ^'Dm‘^m\\LP(n) ^ £0- 


(2.25) 


We extract another subsequence such that W-v^Vm weakly converges to some G in U‘{Q,)'^, and, using the compactness 
of (T’m)m6N and (T>)„)mg]Nr, H-D^Urn — > u in U‘[Q,) and Yi^^Vrn — >■ V* in L^{Q,). Passing to the limit in (2.25) we 
find ||u — i’*||LP(n) > eo- Extending the functions V'D^Vm, and H^^Um by 0 outside fl, we see that, for any 


/r^ 


{yv^Vm{x) ■ (p{x) + n®^Um(a;)divi^(a:)) da: 


and 


/ (Vi7„Um(a:) • P>{x) + ni7,^w„^(a;)divi^(a;)) da: 

jRii 




< W-D^itpln)- 


By limit-conformity of both sequences of gradient discretisations, we can let m —>■ cxo and we find 


IR‘‘ 


( 

Jr'I 


(G • </j(a:)-|-u*(a:)div((5(a:)) da: = / (G • v^(a:)-|-u(a:)div(p(a:)) da: = 0. 


This proves that v,v* G H)|’^(f2) and that G = Vu = Vu*. Poincare’s inequality then gives v = v*, which contradicts 
ll'W — «*IUp(n) > eo. Therefore the sequence (a;m)meN defined by (2.24) satisfies (2.22). □ 


2.5 Polytopal meshes and discrete functional analysis 

Although gradient discretisations are not limited to mesh-based methods (for example it is easy to include spectral 
methods in this framework), a large number of schemes for (2.1) are built on meshes. 

Definition 2.21 (Polytopal mesh): Let fl be a bounded polytopal open subset of R'* (d > 1). A polytopal mesh of O is 
given by T = {M,S,V,V), where: 
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1. Af is a finite family of non empty connected polytopal open disjoint subsets of Q (the cells) such that = UkgmK. 
For any K £ M, \K\ > 0 is the measure of K and hx denotes the diameter of K. 

2. £1 is a finite family of disjoint subsets of (the edges of the mesh in 2D, the faces in 3D), such that any o £ £ is 

a non empty open subset of a hyperplane of R'* and a C fl. We assume that for all if £ Af there exists a subset 
£k of £ such that dK = We then denote by Af^ = {K £ Af : n £ £k}- We then assume that, for all 

(7 £ £, Ma has exactly one element and a C d£l, or Afo- has two elements and a C fl. We let £int be the set of all 
interior faces, i.e. a £ £ such that o d £1, and fext the set of boundary faces, i.e. a £ £ such that a C d£l. For 
o £ £, the (d — l)-dimensional measure of cr is |(t|, the centre of gravity of a is Xa 

3. V = {xk)kgm is a family of points of £l indexed by Af and such that, for all K £ Af, xk € K (xk is sometimes 
called the “centre” of K). We then assume that all cells K £ Ml are strictly a:K-star-shaped, meaning that if x is 
in the closure of K then the line segment [xk, x) is included in K. 

4. V is the set of vertices of the mesh. The vertices that belong to K, for K £ Ml, are gathered in Vk', the set of 
vertices of a £ £1 is denoted by Va- 

For all K £ A4 and for any a £ £k, we denote by the (constant) unit vector normal to <j outward to K. We also 
let dx.o- be the signed orthogonal distance between xk and a (see left part of Fig. 1), that is: 

dx,<T = {x — Xk) ■ n.K,a , 'ix £ a (2.26) 

(note that {x — xk) ■ nx.o- is constant for a; £ cr). The fact that K is strictly star-shaped with respect to xk is equivalent 
to (Ik, a > 0 for all a £ £k- For all TG £ Af and cr £ £k, we denote by Dk,<t the cone with apex xk and base cr, that is 
DK,cr = {tXK + (1 — t)v : £ £ (0,1), y £ a}. The diamond associated to a face cr £ £1 is = UxeA 4 DK,cr- 
The size of the discretisation is Hm ~ sup{/ix : K £ A4} and the regularity factor 6j- is 

Or = max I £ Af , cr £ fxl -h max I ^ : cr £ £!;„* , Ml^ = {K, L}1 . (2.27) 

( dif.o- |-Dx,ct| J [ CtL.CT J 



Fig. 1: A cell AT of a polytopal mesh (left). Two neighbouring generalised hexahedra (right). 


Remark 2.22 (Generalised hexahedra): This definition covers a wide variety of meshes, including those with non-convex 
cells and cells sharing more than one face; in particular, “generalised hexahedra” with non planar faces can be handled; 
such cells have 12 faces (if each non planar face is split in two triangles), but only 6 neighbouring cells. See right of Fig. 
1 . 

Remark 2.23: Since miuo-gfj^ dK,cr is smaller than the radius of the largest ball centred at xk and contained in K, an 
upper bound on 6j- imposes that the interior and exterior diameters of each cell are comparable. 

We now introduce a “polytopal toolbox”, used in the statement of discrete functional analysis results. 

Definition 2.24 (Polytopal toolbox): Let be a bounded polytopal open subset of {d > 1) and let T be a polytopal 
mesh in the sense of Definition 2.21. The quadruplet {X-r,Q, Vr, || • ||t,o.p) is a polytopal toolbox if: 

1. the set Xt-^o is the vector space of degrees of freedom attached to cells and edges (with homogeneous Dirichlet 
boundary conditions): 

X‘K,Q {u {,{Vk) K g , i'^tr) (7(^s) • ^K G R. , V(j £ R. , V a 0ifrr£ £lext} • (2.28) 

2. The mapping IIt- : X'r,o —>■ L°°(£l) is defined by 


Vu £ Xr,o, VA £ Af, Htv = vk on K. 


(2.29) 
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3. The discrete gradient V 7 - : X-j-fi 1 —is defined by 

yK G M, Xrv = 7 ^ ^ kl(w<T - VK)nK,tT = ^ \a\v„nK,<T on K 

(the second equality follows from the property klnx.a = 0, a consequence of Stokes’ formula). 

4. The space is endowed with the following discrete norm, for some p £ (1, 00 ): 

MvgXtp, lkllr.o,p= X] 

K&M CTe£_fc 


Vcr - VK 
dK,f7 


(2.30) 


(2.31) 


In the sequel, T refers to both the polytopal mesh and to the quadruplet (Xr^o, Ht-, Vr, || • ||t,o,p)- 
The discrete gradient Vr satisfies, thanks to Holder’s inequality and to kktc.o- = d\K\, 

\\'^Tv\\LP^Q)d < d p ||u||r,o,p- (2.32) 

Remark 2.25: Note that a polytopal toolbox is not a gradient discretisation, since ||Vr • ||LP(n)‘i i® not a norm on X-r^o 
(consider v € Xr.o such that Va = 0 for all u £ T but vk k 0 for some K G Ml). 

The following lemmas, whose proof can be found in [33, 27], are used to establish Proposition 2.31 below. 

Lemma 2.26 (Discrete Poincare inequality): Let T be a polytopal toolbox of in the sense of Definition 2.24, and let 
9 > Op. There exists C 2 only depending on fl, 9 and p such that for all v G Xpfi we have ||n 7 -u|kp(n) < 6721 |w||r,o,p- 

Lemma 2.27 (Discrete Rellich theorem): Let {Tm)m.eN be a sequence of polytopal toolboxes of D in the sense of Dehnition 
2.24, such that (6'r„)mGN is bounded. If Vm £ Xr„,o is such that (||wm||rm.o,p)meN is bounded, then (nr„Um)meN is 
relatively compact in L^{Q). 

Lemma 2.28 (Discrete approximate Stokes formula): Let T be a polytopal toolbox of L2 in the sense of Definition 2.24. If 
(f G Cr(R.‘^)‘^ and v G Xr.o, then 

/ [Xtv{x) ■ ip{x) + \l'rv{x)AiYtp{x)\ da; 

Jn 

Moreover, if {vcr)crG£K are the exact values at {xcr)cr€£j( of an affine mapping L, then Vpv = VL on K. 

The preceding discrete functional analysis results are useful for the analysis of a wide number of numerical methods, 
thanks to the notion of control of gradient discretisations by polytopal toolboxes. 


< (d|n|) p ||<p|lrvi.oo(Hd)d||u||r.o.phAi- 


(2.33) 


Definition 2.29 (Control of a gradient discretisation by a polytopal toolbox): Let D be a bounded polytopal open subset of 
R'* (d > 1), let ID be a gradient discretisation in the sense of Definition 2.1, and let T be a polytopal toolbox of in the 
sense of Definition 2.24. A control of D> by T is a linear mapping $ : X-d,o —t Xpfi. We then dehne 


l^lk.r = max 


ll^k)llr.o,p 


Ijn®-!; — IIr^(u)||LP(n) 


w (0,r,^)= max 


\\Xt>v\ 


LP(C 2 )‘‘ 


w (D>,T, ^)= max -—-r 

"6^-0,o\{0} ||Vl2U| 


/ [Vi 2 u(a;) — Vr^(n)(®)] da; 
Jk 


LP(n)<^ 


In most of the examples of gradient discretisations in Section 3, the following definition and proposition are used to 
establish the coercivity, compactness and limit-conformity. 


Definition 2.30 (Regularity of a sequence of polytopal meshes): A sequence of polytopal meshes (71n)mgN in the sense of 
Dehnition 2.21 is regular if (0Tm)mGN is bounded and if —>■ 0 as m —>■ -Loo. 



3 Review of gradient discretisations 


12 


Proposition 2.31 (Properties of gradient discretisations controlled by polytopal toolboxes): Let be a bounded polytopal 
open subset of R'* (d > 1). Let (Vm)meN be a sequence of gradient discretisations in the sense of Definition 2.1, and let 
{'Tm)meN be a sequence of polytopal toolboxes of O in the sense of Definition 2.24 such that the corresponding sequence 
of polytopal meshes is regular in the sense of Definition 2.30. We take, for all m € N, a control of by Tm in the 
sense of Definition 2.29, and we assume that 

There exists Cctri > 0 such that, for all m G N, \\^m\\'D^,Tm ^ C'ctri, (2.34) 

lim w"(D™, T™, ^™) = 0, (2.35) 

m—^oo 

lim ^m) = 0. (2.36) 

m—f-oo 

Then (Vm)meN is coercive in the sense of Definition 2.2, limit-conforming in the sense of Definition 2.4, and compact in 
the sense of Definition 2.5. 

Proof We let ujm = max[a;^(r>m, 7)n, ^m), 7)n, ^m)] and M = maxmgNa)m. 

COERCIVITY: using (2.35) and Lemma 2.26, we observe that, for any v € 

l|ni>„'y||LP(n)ti < -^^l|V'D„u||2,p(n)d -P ||nr„$m(u)||iP(n) < AL||Vi?„u||^p(Q)d -P C2||^m(u)||r.o.p- 

Property (2.34) therefore give ||ni 2 „,u||p^p(f 2 )<i < {M + C' 2 C'ctri)| | V-d„u|| iP(f 2 )d, and the coercivity follows. 
Limit-conformity: as stated in Remark 2.8, since is polytopal and therefore locally star-shaped, we only need to 
consider (p G (7(?°(R‘*)'*. By the triangular inequality, (2.35) and (2.33), we have 


J ' ¥’(*) + ^T)rnV{x)A\Ytp{x)^dX 

/ [V'E’„u(a:) — Vrm^ m («)(®)] • ‘Pix)dx 

Jn 


< 


< 


+ l|div<p|| LP' llVr, 


/ [Vr™^m(w)(a;) • (fix) + IlT^^m{v){x)div(fi{x)] da: 
Jn 

/ [^'DmV{x)— Vrm.^ m («)(®)] • <p(*)da; 

Jn 


+ l|div(p|| LP' lIVi, 


+ (d|D|)^ (7,^11$ m («)llr™,o.phA4„ 


(2.37) 


where — 11<(3||^yl,oo(J^d)d. We define tpK = <p(®)da; and notice that \pk\ < C^p and \ip{x) — tpK\ < 

for all X G K. Therefore, since Xj-^^rniv) is constant in each cell. 


[Vi 2 „u(a:) — Vrpi^m(u)(a:)] • (fi{x)dx 


[ [Vi?„u(a:) - Vr™^m(w)(a;)] • (^(a:)d; 

y~^ ( / VT>r,^v{x) ■ [(p{x) - (fiK]dx + <Pk ■ / (VD„u(a:) - Vr„^m(u)(a:))da:) 

KeMrr, ' 

<C^ yy (hMrn. / \^T>mV{x)\dx + I {V (x) — V ^ miv) {x))dx j 

\ Jk Jk j 


K^Mt, 


<C^(yiM^\^\ ’’ + tjJ m ) IIVI2 

We used Holder’s inequality and (2.36) in the last line. Plugged into (2.37) and using (2.34) this gives 


/ p —1 \ p —1 

” -P CUm j -P ||divy 3 ||^p/(f,^Wm-p (d|D|) P C'<^C'ctrlhA 1 „- 

The limit conformity of (Dm)meN follows. 

Compactness: by (2.34), if (V'D„Um)meN is bounded in L^(D)‘^ then \ \^m{vm)\\%r,,o,p is bounded. Applying Lemma 
2.27, we obtain the relative compactness of (nrm$m(u))mgN in L^(D). Since WU-D^Vm — nrm^m(u)||ip(n) —>■ 0 as 
m —>■ 00 by (2.35), we deduce that (nr>,,^Wrn,)mGN is relatively compact in L^{Q). □ 


3 Review of gradient discretisations 

We now study a number of known methods among hnite element, finite volume methods, mimetic methods and related 
discretisation schemes which are all based on polytopal meshes. Each of the following sections is devoted to a particular 
method which is shown to be the gradient scheme of a gradient discretisation referred to as D; for each method we define 
a regular sequence (I’m)meN of gradient discretisations, based on the method itself and on the regularity of a polytopal 
mesh (Definition 2.30), and we show the following property. 
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Fig. 2: Partitions for mass-lumping of the Pi (left) and non-conforming Pi (right) finite element methods. 


The regular sequence (22m)meN is coercive, consistent, limit-conforming and compact in the sense 
of the definitions in Section 2.1. 

The proof of {V) relies on the notions of LLE gradient discretisations (section 2.2), barycentric condensation (section 
2.3), mass lumping (section 2.4) and polytopal toolbox (section 2.5). 


3.1 Pfc finite element methods 

3.1.1 Conforming methods: Pfe finite elements 

Let T be a simplicial discretisation of fi, that is a polytopal discretisation in the sense of Definition 2.21 such that for any 
K € M we have C&t:A{Ek) ~ d-|-l. Let k £ N\{0}. We follow Definition 2.11 for the construction oiV — {X-pfi, Vn, IIxi) 
by describing the partition of Q, the functions ai and the local linearly exact gradients in the elements of the partition. 

1. The set I of geometrical entities attached to the dof is 7 = and the set of approximation points is S = I, 

where = Uksm ^ points x of the form 

x= with (i„)„gvji: €{0, such that = k. (3.1) 

Then In = (the subset of the interior vertices) and Ign = (boundary vertices), and the partition of Q is 
given by U^M. For f/ = 7L G W, we let lu = 

2. The reconstruction operator Dn in (2.13) is defined using the basis functions (Q;v)^gy(fc), called in this particular 
case the Lagrange interpolation operators and defined the following way: in each cell K, is the polynomial 
function of x with degree k, such that av(v) = 1 and av(v') = 0 for all v' G \ {v}. This leads to 

Vn G X-pfi, \/x G n, Ti'Dv{x) = ^ Vya^{x). 

3. The linearly exact gradient reconstruction in K is 

'ix G K, Gkv(x) = UvVav(a:) = V(ni3ii)(a:). 

vev)C) 

4. We have Yl-pv G ITq’^)!!) so the Poincare inequality in lL))’*’(fl) implies that ||Vc • \\LP(n)d is a norm on Xt>,o- 

Defining the regularity of a sequence of P* discretisations {T>m)m.eN merely as the regularity of the underlying polytopal 
meshes (Definition 2.30) is sufficient to obtain the boundedness of (regLLj;(Ilm))meN; hence Proposition 2.14 implies the 
consistency. Since XH-dv = Xt>v we have Wt> = 0 and the limit conformity is trivial; the coercivity and the compactness 
are consequences of the Poincare inequality and the Rellich theorem in respectively. This establishes {V) for 

Pfe gradient discretisations. 

3.1.2 Mass-lumped Pi finite elements 

We construct a mass-lumped version of the Pi gradient discretisation as per Definition 2.19, with the natural geometrical 
entities attached to the dof = V. Subdomains (K)vev with points of V as centres can be constructed in various 
ways. One way is to define K, as the set of all y G 17 such that ay(y) > ayi{y) for any other v' G V. The left part of 
Figure 2 illustrates the construction of the partitions (K)vev in the case d — 2 (then (K)v6V is sometimes called the 
barycentric dual mesh of T). 

A Taylor expansion in each K H 77 then shows that Estimate (2.22) holds with u)rn = hM-m, and thus by Theorem 2.20 
we see that the mass-lumped Pi gradient discretisation satisfies the property {V). 
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3.2 Non-conforming Pi finite elements 

3.2.1 Standard non-conforming Pi reconstruction 

Non-conforming Pi finite elements consist in approximating the solntion to (2.2) by functions that are piecewise linear on 
triangles and continuous at the edge midpoints - but not necessarily continuous on the whole edge. These approximating 
functions therefore do not lie in Hq{Q,), and do not satisfy the exact Stokes formula; hence the name “non-conforming”. 
Let T be a simplicial mesh of that is a polytopal mesh in the sense of Definition 2.21 such that for any Lf G Af we 
have Card(T_K') = d+ 1. We refer to Definition 2.11 for the construction of D. 

1. The set of geometrical entities attached to the dof is I = £ and the approximation points are S — {xa)cr€S- Then 
Ici ~ Tint and Ion = Text, and the partition of Q is given by M = M. For all [/ = K gU, we let lu = £k- 

2. The reconstruction IIt> in (2.13) is defined using the affine non-conforming finite element basis functions {acr)ae£ 

defined by: is linear in each simplex, ao-(®o-) = 1 and aa{x„i) = 0 for all a' G T\{(t}. This leads to 

Vu G Xt>,o, Va: G fl, XIt>v{x) — ^ Uo-Qo-)®). 

(7^S 

3. The linearly exact gradient reconstruction in K is defined by the constant value 

yx € K, Qkv(x) = ^ Va\7a„-(x). 

4. The fact that IjVrj • ||i,p(n)d is a norm on Xti,q is deduced from the injectivity of the mapping defined in the 
course of the proof of the property {V) below. 

The regularity of the non-conforming Pi gradient discretisations is then defined as the regularity of the underlying 
polytopal discretisations (7)n,)meN (see Definition 2.30). 

Proof [Proof of the property (V) for non-conforming Pi gradient discretisations] 

We drop the index m from time to time for sake of legibility, and all constants below do not depend on m or the 
considered cells/edges. Let us define a control of D by T in the sense of Definition 2.29, where T is the simplicial mesh 
associated to D, with xk ~ xk = ^o-gSk centres of gravity of the cells K. We define the linear (injective) 

mappings $ : Xt>^,o — Xr^^ by ^{u)k = ^ ^ ^vu^xk) and ^(u)^ = = IIt>u{x^). 

Since ^{u)k = I1t>u{xk) and Qku = V(niiri) in K, we get 


$(u)ct - ^{u)k = Gku ■ {x„ - Xk)- 


Therefore, since < 6r, 

’ d.K,„ — dK,„ — ’ ’ 


^ |cr|dx,a 
ct^Sk 


^(m),^ - ^{u)k 

dK,a 


< e^^d\K\ \gKu\^. 


(3.2) 


This implies (2.34). We now observe that the affine function reaches its extremal values at the vertices of K. It is 
easy to see that acr{ycr) = 1 — d, where Vo- is the vertex opposite to the face o, and that q;ct(Vo.') = 1 for all a' ^ a. 
Therefore, for x G K, 


\IItju{x) - ^(u)k| 


E 

CT^Sk 


^{u)K)a-a{x) 


< (d -I- 1) max(l, d — 1) max \Qku ■ (xa — xk)\- 
creSK 


This inequality implies cu'^(iD,T, ^) < {d + l)max(l,d— l)hM and therefore (2.35) holds. Finally, recalling that IIt)U 
is affine in each simplex K and that Vr is exact on interpolants of affine functions (cf. Lemma 2.28), we see that 
V-dw = Vr^(w) in £l. Hence tu^(iD,T, ^) = 0 and (2.36) holds. Proposition 2.31 therefore shows that ('Dm)meN 
is coercive in the sense of Definition 2.2, limit-conforming in the sense of Definition 2.4, and compact in the sense of 
Definition 2.5. 

Since non-conforming Pi gradient discretisations are LLE gradient discretisations, the consistency of (Dm)mGN follows 
from Proposition 2.14 by noticing that is controlled by □ 


3.2.2 Mass-lumped non-conforming Pi reconstruction 

Let us recall that, in the case d = 2, for all pair (a, a') G T^ with a yf a', there holds 

/ aa-(x)a„.>(x)dx = 0. 

Jn 

This property ensures that the non-conforming Pi method has a diagonal mass matrix, which is useful for computing 
(2.10). However, Property (2.9) is not satisfied, which might prevent the usage of the non-conforming Pi scheme for some 
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nonlinear problems. To recover a piecewise constant reconstruction, and thus (2.9), we apply to the preceding gradient 
discretisation the mass lumping process as in Definition 2.19. Recalling that the set of geometrical entities attached to 
the dof is I — we define the subdomains (Vi)ig/ of Definition 2.19 as the diamonds {Dcr)cr€e around the edges. The 
right part of Figure 2 illustrates the construction of this partition. Since IlxiV is linear and V(ni 5 u) = V-pV in each cell, 
and since = IlTiv{xcr) on D^, an order one Taylor expansion immediately provides Estimate (2.22). Property {V) 

for the mass-lumped non-conforming Pi gradient discretisation is then a consequence of Theorem 2.20. 


3.3 Mixed finite element MT^ schemes 


The RTfc method is the only one presented here for which the gradient discretisation is only constructed for p — 2. All 
the other gradient discretisations are constructed for any p G ( 1 , oo). 

Let T be a simplicial discretisation of D as for the non-conforming Pi scheme. We fix fc G N and introduce the following 
spaces 

Vh = {we iL\n)f : w\k G KTfe(A'), VfF G M}, Vf" = V;, n JTdiv(D), 

Wh = {peL\n) : GPfe(/r), VKgM}, M^= i/i : U GPfc(f"),f‘|an = 0 

I aGS 


where 

• Pfc(A') is the space of polynomials of d variables on K of degree less than or equal to k. 

• Pfc(cr) is the space of polynomials of d — 1 variables on a of degree less than or equal to k. 

• RTfc(A') = Pfe(A')‘* + x¥k{K) is the Raviart-Thomas space of order k defined on K. Here, Pfe(A') C ¥k{K) is the 
set of homogeneous polynomials of degree k. 

We construct a gradient discretisation (for p = 2 only) inspired by the dual mixed finite element formulation of Problem 
(2.1) as in [8]. Assuming that A is constant in each cell K, the dual mixed finite element formulation of (2.1) is 

iv,q,X) eVhXWhX Ml 

f w{x)-A{x)-^v{x)dx- f g(a:)divm(a:)da:-b ^ [ \{x)w\k{x) ■ nK,crds{x) = 0, Xw e Vh, 

J K J K J(7 


L 


7^Sk ' 

ip{x)divv{x)dx = / 'ip{x)f{x)dx, Vi/> G Wh,XK G M, 


p{x)v\k{x) ■ nK,ads{x) + / pl{x)v\l{x) ■ nr.,ads)®) = 0, Vcr G Tint with Ma = {K, L}, V/r G M°. 
J <y 


(3.3) 


We again refer to Definition 2.11 for the construction of D = (Xri^o, V-d, Hi:)). We consider the standard basis 

of Wh, and (G')jer“ the standard basis of Ml These two standard bases are respectively associated to the set of points 
located in the cells and the set of points located on the faces of the cells. These points are defined in a similar 
way as (3.1). 

1. The set of geometrical entities attached to the dof is 7 = U 7^, and the set of approximation points is also 
S = 7'^ U 7". Then 7n = 7^^ U 7^^ and Ign = where 7;" = 7" n H and 7e" = 7" n 90. The partition of O 
is given by 77 — M. We denote by ll the set of all points of which are in K, and by the set of all points 
of which are in a. Then, for all 7/ = 77 G 77, 7;/ = U Uct 6 £^ ■ 

2. The reconstruction (2.13) is applied with a; = ■i/ji for all i G ll and at = 0 for all i G Uo-efj^ ■ This leads to 

Vu G X-pfi, Xx G O, II'dv(x) = ^ Vi'tpi{x). 

3. For all K G M, the linearly exact gradient reconstruction is locally defined in K by: Qkv is the function such that 
AQkv G RTfc(77) and 


Viipi{x) I diviy(a 3 )da; 

/ 

\jelM J 

4. If II VculliPfnjd = 0 then Gku = 0, and {v,q,X) defined by v\k = AQku, q = J2ieiw Uiipi, A = ^ 

solution to (3.3) with f — 0. The invertibility of this system implies that q — 0 and A = 0, and therefore Ui — 0 
for i G 7'^ and Uj = 0 for j G 7^. Therefore |j V-d • || LP(n)d' is a norm on X-d,o- 
The proof of the equivalence between the corresponding gradient scheme (2.3) and the Arnold-Brezzi mixed hybrid 
formulation (3.3) is found in [34], along with the proof of the property (V) for a regular sequence of polytopal meshes in 
the sense of Definition 2.30. Note that, in the case k — 0, the property of piecewise constant reconstruction holds. 


ViyGRTfc(77), / w(x) ■ QKv(x)dx + 


Jk 
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3.4 Multi-point flux approximation MPFA-0 scheme 

We consider in this section two particular cases of the MPFA-O scheme [1]. They are based on particular polytopal 
meshes of Q in the sense of Definition 2.21: Cartesian for the hrst case, and simplicial for the second case. In each of 
these cases, for K G A4 we let xk = xk be the centre of gravity of K and we dehne a partition {VK,v)\ieVK ^ 
following way (see Figure 3): 

• Cartesian meshes: Vk,v is the parallelepipedic polyhedron whose faces are parallel to the faces of K and that has 
xk and v as vertices. We define, for a G £ and v G Vo-, ®(o,v) = ^o (note that these points are identical for all 
V G Vo, see Remark 2 . 12 ). 

• Simplieial meshes: we denote by {x))-^^\>k barycentric coordinates of re in if (that is x—xk = 

Xk), (x) > 0 and P^(^) ~ 1) ''^® define Vk,v as the set of a; G if whose barycentric coordinates 

(/3,(f(a:))„/gVjf satisfy (x) > P^ (x) for all v' G Vk \ {v}. For a G £ and v G Vo, ®(o,v) is the point of a whose 
barycentric coordinates in a are (a:(o,v)) = l/(d+ 1) for all v' G Vo \ {v}, and /3v (a;(o,v)) = 2 /(ci+ 1 ). 




Fig- 3: Notations for MPFA-0 schemes defined on Cartesian (left) and simplicial (right) meshes. 


We then follow the notations in Definition 2.11 to construct the MPFA-0 gradient discretisations in both cases: 

1. The set of geometrical entities attached to the dof is 7 = At U {(cr, v) : a G £, v G Vo} and the family of 

approximation points is S' = ((a:x)if6Ai, (a:(o,v))oG£,veVo)- We define 7n = Af U {(cr, v) : a G Sint, v G Vo} 
and Jen = {(c,v) : a G Sext, v G Vo}. The partition is 77 = (VK,v)KeM,veVK ■ For any U = Vk,m, we set 

£k,v = {<j G £k : V G Vo} and lu = {K} U {(cr, v) : a G S/c.v}. 

2. The functions Oi are dehned by at = 1 for i = K and Oi = 0 for i = {a, v), which means that 

Vu G X-Dft , V7f G AI, Va: G 7f, Ii'Dv{x) = vk- 


3. Setting a„ = Vk,v H cr, the gradient reconstruction on [/ = Vk,v is 


Va: G Vk,v , Gvk,„^{x) 


1 


|crv|(w(o,v) - nK)nif.o. 


4. As in the case of the non-conforming Pi element, the fact that HV-p • ||iP(n)<i is a norm on Xxip is deduced from 
the injectivity of the mapping $ defined in the course of the proof of the property {V) below. 

For such a gradient discretisation, the gradient scheme (2.3) is a finite volume scheme. Indeed, by selecting a test function 
with only non-zero value vk = 1 in (2.3), we obtain the flux balance 

E E FK,<T,^iu) = / f{x)dx, where FK,a,y{u) = / Qvk,„'u{x) ■ nK,ad'y{x). (3.4) 

o-e£jfvev„ “'re “'o-„ 


Selecting a test function with only non-zero value U(o-,v) = 1 in (2.3) leads to the conservativity of the fluxes: 

FK,a,y(u) -I- FL,cr,y{u) = 0 for all cr G Sint with Mcr = {K, L}, and all v G Vo-. (3.5) 


We can also locally express the degree of freedom U(ct,v) in terms of (uk) ifiveVj^ • For a given v G V this is done by 
solving the local linear system issued from (3.5) written for all cr such that v G Vo- After these local eliminations of 
U(o,v), the resulting linear system only involves the cell unknowns. This discretisation of (2.1) by writing the balance 
and conservativity of half-fluxes T^/c.o.v, constructed via a local linearly exact gradients, is identical to the construction 
of the MPFA-0 method in [1]. This demonstrates that the gradient discretisation constructed above indeed gives the 
MPFA-0 method when used in the gradient scheme (2.3). 
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Remark 3.1: The identification of MPFA-O schemes as gradient schemes is, to our knowledge, restricted to the two cases 
considered in this section (Cartesian and simplicial meshes). In the case of more general meshes for the approximation of 
(2.1), the discrete gradient defined by the MPFA-O scheme can be used in the finite volnme scheme (3.4)-(3.5); however, 
the gradient scheme (2.3) built upon this discrete gradient cannot be expected to converge, since the corresponding 
gradient discretisation may fail to be limit-conforming and coercive. 

The regularity of the MPFA-O gradient discretisations is defined as the regularity of T (see Definition 2.30). Although 
references [41, 40] include proofs of the property ("P), let us show how the generic tools presented in Section 2.5 enable 
very quick proofs of this result. 

Proof [Proof of the property {V) for MPFA-O gradient discretisations] 

We drop the indices m for sake of legibility. We consider the polytopal mesh T = (Af,T^'P, V^) where the sets (M,P) 
are those of the original polytopal mesh, £' = {uv ; cr G T, v G V}, and V' is the set of all vertices of the elements of 
. We define a control of D by T (in the sense of Definition 2.29) as the isomorphism $ : Xxi,q — X-r,o given by 
^{u)k = uk and $(m)ct„ = «(<,.,v)- We observe that 



with Cs = 1 for parallelepipedic meshes, and Os > 0 depends on an upper bound of the regularity of the mesh for 
simplicial meshes. Therefore II > C'3||^('u)|lr,o,p (2.34) is proved. Since II-dm = nr^(M), we get 

0 J^('D,T, $) = 0, which proves (2.35). Finally, we have 


/ V-D«(a;)da: - E E = |a'|($(u),,-$(u)K)nxx = |If| Vr^(n)|K. 

^ o’Cfp'veVcr 


This shows that = 0, which establishes (2.36). Proposition 2.31 therefore shows that (iDm)meN is coercive 

in the sense of Definition 2.2, limit-conforming in the sense of Definition 2.4, and compact in the sense of Definition 2.5. 
It is proved in [41, 40[ that the definitions of the approximation points S give the LLE property in both the Cartesian 
and simplicial cases. Hence, the consistency of (iDm)m6N follows from Proposition 2.14. □ 


3.5 Discrete duality finite volumes 

The principle of discrete duality finite volume (DDFV) schemes [44, 22, 12, 45, 4, 17[ is to design discrete divergence 
and gradient operators that are linked in duality through a discrete Stokes formula. Since discrete operators and an 
asymptotic Stokes formula are at the core of gradient schemes (see Definitions 2.1 and 2.4), it is not a surprise that they 
should contain DDFV methods. This was already noticed without proof in [36]; we give here a precise construction and 
proof of this result in the 3D case. Note that the same tools can be used for the 2D case [27]. 

Two 3D DDFV versions have been developed: the CeVe-DDFV, which uses cell and vertex unknowns [46, 20, 5], and 
the CeVeFE-DDFV, which uses cell, vertex, faces and edges unknowns [17, 18[. The coercivity properties of the two 
methods differ: the CeVe-DDFV does not seem to be unconditionally coercive on generic meshes, whereas the CeVeEE- 
DDFV scheme is unconditionally coercive [23]. We show here that this latter method is a gradient scheme. To do 
so, we introduce a gradient discretisation on a general octahedral mesh (possibly including degenerate cells), and we 
show that when the octahedral cells of this mesh are the “diamond cells” of a CeVeFE-DDFV method, the gradient 
scheme corresponding to this gradient discretisation is the CeVeEE-DDFV scheme. The standard CeVeFE-DDEV scheme 
corresponds to hexahedral meshes, which can be seen as degenerate octahedral meshes (each cell has six vertices, but 
three of them are aligned so only six physical faces are apparent). Although it was known to specialists that, as done here, 
the construction could be performed on general octahedral meshes (and that the corresponding DDEV method satisfies 
the discrete duality formula [7]), this was not reported before. It should also be noticed that our presentation gives a 
complete description of the CeVeFE-DDFV method using only one mesh instead of the usual four meshes. As shown in 
[6, Section IX.B] for 2D DDFV methods, the other three meshes can be reconstructed from the ocatahedral (“diamond”) 
mesh. However these three meshes are not used to construct the method here. The vision of DDFV methods based 
solely on one mesh (the “diamond” mesh, or octahedral mesh here) actually corresponds to the vision adopted in the 
implementation of the schemes. 

Let T = {M,£,P, V) be a polytopal mesh of D in the sense of Definition 2.21, such that the elements of A4 are octahedra 
(open polyhedra with eight triangular faces and six vertices, not necessarily convex; five vertices may be coplanar), and 
the element of £ are the triangular faces of the elements of M. Each £k has 8 elements, each Vk has 6 elements, and 
each Vct has 3 elements. For any K G M, the centre of K is defined by xk = \ X^vgVk Definition 2.11 to 

construct an “octahedral” gradient discretisation T> = (Vu^o, V-d, H-d) (see Figure 4 for some notations): 

1. The set of geometrical entities attached to the dof is / = V, the set of approximation points is S' = V. We let 
7n = V n D and Iga = V n dQ,, and the partition is W — M. For U = K £14 we define lu = Vk- 
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2 . For K £ M and v £ Vk, we denote by Vk,m the octahedron formed by xk, v, and the four other vertices of K 
that share a face of K with v. We then use the reconstruction (2.13) with the functions (av)verjc defined by 

Va; £ K, Wv £Vk , a^{x) = ^XVk.v(®). (3-6) 

where Xv^ v denotes the characteristic function of Vk,v This leads to 

Vu £ Xt,,o ,yK£M,yx£K, nT,u{x) = (3.7) 

veVjc 


3. For K £ M and u £ Xufi, the cell gradient is defined by 

Va; £ K, Qku{x) = ^ \G\u„\\K,a , where Mo- = ^ ^ for all a £ £k- (3.8) 

vgVct 

4. The proof that ||Vi5 • ||LP(n)rf is a norm on Xt>,o is done in Lemma 3.4. 




Fig. 4: Left: octahedral cell K for the CeVeFE-DDFV scheme. Center: illustration of Vx^y Right: construction 
of a degenerate octahedron from a non-conforming hexahedral mesh in a heterogeneous medium {CDEF 
is the intersection of the boundaries of two non-matching hexahedral cells). 


Remark 3.2 (Octohedra and heterogeneous media): An octahedral mesh of can be obtained in the case where the domain 
XI is the disjoint union of star-shaped hexahedra, considering the octahedra obtained from the centres of neighbourhing 
hexahedral cells, and the four vertices of their interface. This also works for non-conforming hexahedral meshes, for 
which interfaces between cells may be different from the physical faces of the cells. 

In the case of a heterogeneous media, in which the material properties (e.g. the permeability A in (2.1)) are constant 
inside each hexahedral cell but may be discontinuous from one cell to the other, it is usually preferable to construct 
octahedral cells that are compatible with these heterogeneities (i.e. such that the material properties are constant inside 
each octahedron). This prevents from introducing a non-physical average of A in the gradient scheme (2.3), which would 
lead to a loss of accuracy of the approximate solutions. Such an octahedral mesh can be constructed fairly easily as 
illustrated in Figure 4 (right). Each of these octahedra is built from the centre of an hexahedral cell, the four vertices of 
the interface between this cell and a neighbouring hexahedral cell, and a point selected on this interface. This interface 
need not be planar. 

Remark 3.3 (Reconstruction operator): The choice (3.6) of n® is driven by our desire to construct a gradient discretisation 
whose gradient scheme is exactly the CeVeFE-DDFV method, for particular octahedral meshes. This choice ensures that 
the discrete duality formula holds true but, as explained in Section 2.1 (and as already noticed in [3, Appendix C] for 
CeVe-DDEV methods), it is not adapted to certain non-linear models. 

The following lemma proves that the previous construction gives an LLE gradient discretisation. It will also prove useful 
to show that this gradient discretisation gives back the CeVeFE-DDEV method, and to establish Property {V). 

Lemma 3.4: Let T> be the octahedral gradient discretisation defined as above. Eor any u £ X-Dfi and any K £ Mi, the 
constant discrete gradient Qku is characterised by 

For all opposite vertices (v, v') of K, Qku ■ (v — v^) = (3.9) 

As a consequence, ||Vi3 . ||i,p(f 2 )d is a norm on X-pfi- 
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Remark 3.5: The opposite vertices in the octahedra in Figure 4 are {A, B), (C, D) and [E, F). 

Proof We first note that, since the three directions defined by the three pairs of opposite vertices in K are linearly 
independent, (3.9) indeed characterises one and only one vector Qku G R®. We therefore just have to show that the 
gradient defined by (3.8) satisfies (3.9). We have 

^ (3.10) 

' ' veVjf aeSKWev^ 

Let us consider for example the case where v = A in Figure 4. For a triangular face a we can write Icrlnic^o- as the 
exterior product of two of the edges of cr (with proper orientation). This gives 

|a|nK.<T = ^(^ X A^ + X Az5 +Az5 X A^ +X AC^) 

o-G£k |vGVcr 

= xaP + ^x A^) = X 

Applying this to all vertices of K, and since 1^1 = |Aif with Ak = det(A^, Cl^, E^), we deduce from (3.10) that 

Gku = ^— (j^^B — ua)c£> X E^ + {ud — uc)E^ x + (mf — ue)A^ x . 

Property (3.9) is then straightforward. Considering for example the case (v, v') = {B,A), the formula follows from 
(I^ X A^) • A^ = (A^ X C]5) • A^ = 0 and (^ x ■ A^ = det(^, W, A^) = Ak- 

Let us now prove that |j Vi> • ||FP(n)£i is a norm on Xt),o- Assume that \\X■du'^lp(ci)’>- ~ 0- Then for any K G At, Qku = 0. 
Let us take a boundary octahedron K. One of its faces a is entirely contained in dO, and thus by the boundary conditions 
all (uv)vev^ vanish. Using (3.9) we infer that all values of u at vertices opposite to Vo- also vanish. Since any vertex in 
K either belong to Vo- or is opposite to a vertex in Vo, this shows that u vanishes on all vertices of K. We then conclude 
by induction on the number of octahedra in the mesh that u = 0. □ 

It is now easy to detail the relationship between the octahedral gradient discretisation T> and the CeVeFE-DDFV scheme. 

Lemma 3.6 (CeVeFE-DDFV is a gradient scheme): For any polytopal mesh T of Q., there exists an octahedral mesh T of 
Q. such that, if T> is the octahedral gradient discretisation defined as above from T, then the gradient scheme (2.3) for T> 
is the CeVeFE-DDFV method on T. 

Proof The CeVeFE-DDFV method on T has cell, vertices, edge and face unknowns, and its discrete gradient is piecewise 
constant on so-called “diamond cells”. A diamond cell is an octahedra as in Figure 4 (left), but with E chosen on the 
segment [A, B] - hence, the octahedra actually degenerates into an hexahedra. The segment [A, B] corresponds to an 
edge of the primal mesh of the CeVeFE-DDFV method, E is a point on a face of this mesh that contains [A, B], and D 
and C are points inside the cells on each side of this face. Let us take T the polytopal mesh of D made of the diamond 
cells (degenerate octahedra). It is proved in [17, Lemma 3.1] that the CeVeFE-DDFV discrete gradient satisfies (3.9); 
hence, this gradient is Vr>. It is then just a matter of applying the discrete duality formula [17, Theorem 4.1] on the 
formulation [17, Eq. (5.4)] of the scheme to see that the CeVeFE-DDFV scheme for (2.1) is indeed the gradient scheme 
(2.3) for V. □ 

Remark 3.7 (CeVeFE-DDFV and heterogeneities): Except on the boundary of the domain, the diamond cells of the CeVeFE- 
DDFV method are always spread on two neighbouring cells. If the primal mesh is aligned with heterogeneities of the 
medium, then these heterogeneities are actually averaged in the formulation of the CeVeFE-DDFV scheme (see [17, Eq. 
(5.2)]). An option to better deal with heterogeneities is to use the m-DDFV method [12, 19]. Additional unknowns are 
introduced and eliminated by writing local flux conservativity. These eliminations are easy to perform for linear models, 
but require solving local nonlinear equations in the case of nonlinear models. 

As shown in Remark 3.2, even starting from a primal hexahedral mesh aligned with the ground properties it is possible 
to construct a (degenerate) octahedral mesh that is also aligned with the heterogeneities of the medium. Using this 
octahedral mesh ensures a better accuracy than the “standard” CeVeFE-DDFV method, at a reduced computational 
cost compared to the m-DDFV method in the case of nonlinear models (no nonlinear local equation to solve). 

Let us now turn to the analysis of the properties of the octahedral gradient discretisation. The regularity of a sequence 
of octahedral gradient discretisations (Pm)meN is merely defined as the regularity of the underlying octahedral meshes 
(7jn.)mGN of D in the sense of Definition 2.30. With this definition of regularity, the property (V) holds for octahedral 
gradient discretisations. 
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Proof [Proof of the property (P) for octahedral gradient discretisations] As usual, we sometimes drop the indices m for 
sake of legibility. Let us define a control : A-d,o —t Xrfi of P by T in the sense of Definition 2.29, where T is the 
octahedral mesh with, for all K £ A4, the centre of K defined by xk ~ \ SveVjf 

VLf £ M , ^{u)k = ^ ^ Uv and Va £ £, ^{u)a = ^ ^ Mv. (3.11) 

v^V<7 

Let a £ £k, and let us denote by (vi)i=i,2.3 the three vertices of o and by v' the vertex opposite to Vi in K. We then 
have, for u £ X-vfi, 

^{u)a — ^(u)k = + Mv 2 + Mva) — ^ (Mvi + Mva + Mvj + Mvi; + Mv' + Mv') 

O 0 12 6 

1111 . ^ ^ 

= g(Wvi - Mv'i) + g(Mv2 - Mv^) + g(Wv3 - Mv^) = q'^QkU ■ (Vi - v'), 

i=l 

thanks to Lemma 3.4. Therefore, since < 6-r, we have X)o-e£if | | < 

and (2.34) follows. 

For K £ M and a.e. x £ K there are three vertices (vi)i=i_2,3 C Vk such that x £ Vk,^^- Denoting by v' the vertex in 
K opposite to Vi, Lemma 3.4 again gives 

1111 . ^ ^ 
ni,n(a;) - nr^(M)(a:) = -(Wvi - Mv'^) + g(«v2 - K 0 + g(«v3 - K'J = (vi - v'). 

i = l 

This shows that |n-Du(a;) —ni2/$(M)(a;)| < ^hM\GKu\ and thus that Ijllort —nr^(M)||LP(n) < I^ArjlV-DMjj^pjnjd. Hence 
CJ^(P, T, $) < and (2.35) holds. The definition (3.8) implies that = V7-€>(u), and therefore u)^{T>, T, $) = 0, 
which implies (2.36). By Proposition 2.31 we deduce that (Pm)mGN is coercive in the sense of Definition 2.2, limit- 
conforming in the sense of Definition 2.4, and compact in the sense of Definition 2.5. 

The bound on 9j- forces the three vectors (v — v')(„ „/) opposite in K lo be of similar length and “really non-coplanar” - 
that is to say with a determinant of order the cube of their similar length - which gives an estimate on ||f/if||tx>. Hence, 
the regularity factor regLLj;(Pm) is bounded and Lemma 3.4 proves the consistency. □ 

3.6 Hybrid mimetic mixed schemes 

3.6.1 Fully hybrid scheme 

Since the 50’s, several schemes have been developed with the objective to satisfy some form of calculus formula at the 
discrete level. These schemes are called mimetic finite difference (MFD) or compatible discrete operator (CDO) schemes. 
Contrary to DDFV methods that design a discrete operators and duality products to satisfy fully discrete calculus 
formula, MFD/CDO methods design discrete operators that satisfy a Stokes formula that involves both continuous and 
discrete functions. Depending on the choice of the location of the main geometrical entities attached to the dof (faces or 
vertices), two different MFD/CDO families exist. We refer to [47] for a review on MFD methods, and to [11, 10] (and 
reference therein) for CDO methods. 

A first MFD method, hereafter called hybrid MFD (hMFD), is designed by using the fluxes through the mesh faces as 
initial unknowns. This requires to recast (2.1) in a mixed form, i.e. to write q — AVu and div(q) = /, and to discretise 
this set of two equations. The resulting scheme takes a form that is apparently far from the gradient scheme (2.3). 
It was however proved in [28] that this hMFD can be actually embedded in a slightly larger family that also contains 
hybrid finite volume (HFV) methods [33] and mixed finite volume (MFV) methods [24]. This family has been called 
hybrid mimetic mixed (HMM) schemes; each scheme in this family can be written in three different ways, depending on 
the considered approach (hMFD, HFV or MFV). The HFV formulation of an HMM scheme is very close to the weak 
formulation (2.2) of the elliptic PDF; it actually consists in writing this weak formulation with a discrete gradient and 
a stabilisation term (bilinear form on (u, v)). It was proved in [29] that the discrete gradient can be modified to include 
the stabilisation terms, and thus that all HMM methods - which means all hMFD methods also - are actually gradient 
schemes. 

The discrete elements that define an HMM gradient discretisation are the following. We again refer to Definition 2.11 
for the construction of V. 

1. Let T be a polytopal mesh of fl as in Definition 2.21. The geometrical entities attached to the dof are I = M U £ 
and the approximation points are S = ((xK)KeM, (xcr)a-e£). We let In — A4 U fint and Ign ~ £ext, and we have 
Xt>^o = as defined by (2.28). The partition is W = {DK,cr : K £ M,a £ £k} and, for U = DK,a £. U, we set 
Iu = {K}\j£k. 

2. The reconstruction (2.13) is defined by the functions Uk = 1 in K and uk = 0 outside K, for all K £ A4, and 
Ocr = 0 for a £ £k- Recalling the definition (2.29) we therefore have 

Vw £ Xt>,o, Vif £ A4, Va; £ K, 11tiv{x) = nr'y(a;) = vk- 



3 Review of gradient discretisations 


21 


3. 


Recalling the definition of the polytopal gradient (2.30), we start the construction of the discrete gradient by 
setting 


Vn G Xt>,o , VR" G A4 , Xkv = (V7'u)|x 


1 

W\ 


|(T|nCTnic,CT. 


(3.12) 


This gradient is linearly exact, but it does not lead to the norm property: it can vanish everywhere even for n yf 0 
(it suffices to take Va = Q for all a G £)■ We therefore add a stabilisation that is constant in each half-diamond, 
and that vanishes on interpolations of affine functions: 


\/ff 

Vu G X-Dfi , VUic.o- G W , Va; G Dk,<t , Qd^ ~ ^kv + - - [£KRK{QK{v))]a-nK,tT, 

' ClK,cr 

where 

• Qk(v) = {V„ - VK)a€SK^ 

• Rk : is the linear mapping defined by Rk{0 = iRK,a-{^))a-e£K with 


(3.13) 




• Ck is an isomorphism of the vector space Im(RK) C R^^. 

4. The norm property of IjV-p • ||LP(n)‘* on X-pfi follows from the estimate (see [29, Lemma 5.3]) 

Ijwjjr.o.p < C'4||Vi7M||p^p(n)d, 


(3.14) 


where C4 > 0 depends on an upper bound on the regularity factor 9j- of the polytopal mesh and on the regularity 
factor ((r>. This last factor is defined as the smallest number such that, for all if G Af and all ^ G R^^, 


(tGSk 


RK,a{0 


d 


•K,a 


< E 1^^ 

(tGSk 


[CkRk{0]c 


d. 


K,a 


< Co E \^K,a 

(tGSk 


Rif.CT(C) 


dp 


Remark 3.8: The face degree of freedom Va- corresponds to the hybridisation of the hMFD methods. 


The freedom of choice of the isomorphisms {Ck)k£M ensures that all hMFD, HFV and MFV schemes are covered by 
the framework (there are several such schemes, due to different possible stabilisation parameters). More precisely, [29] 
proves that for any HMM scheme 5 on T there exists a family of isomorphism (Tx)ifGAi such that, if V is defined as 
above, then the gradient scheme (2.3) is S. 

The proof of the property (V) for HMM methods was originally given in [29]. We show here how the notion of gradient 
discretisations controlled by polytopal toolboxes notably simplifies this proof. We say that a sequence of HMM gradient 
discretisations (iDm)meN is regular if the sequence of underlying polytopal meshes (7)n)meN is regular in the sense of 
Definition 2.30, and if (C'D„)meN is bounded. 

Proof ]Proof of the property (V) for HMM gradient discretisations] Let ('Dm)me'N be HMM gradient discretisations 
built on polytopal meshes (Tm)meKi, and let us define a control of Pm by Tin in the sense of Definition 2.29. We drop 
the index m from time to time. Since X-p^o = X^-^, we can take = Id. Estimate (2.34) is given by (3.14). Relation 
(2.35) follows immediately since w^(P,T, €•) = 0, owing to H-du — Hy-u = Hy-^)!!). Recalling that ]Dk.ct] = we 

have 

[ V-DU{x)dx ^ \K\\/kU+ E W\[^KRK{QKiu))]anK,a- (3.15) 

The definition of Rk and the property \cr\nK,a-{xa- — xk)^ ~ l^jfd (a consequence of Stokes’ formula) show that 

for any rj G Im(Rir) we have ]cr]77CTnx,o- = 0. Hence, since Im(Tx) = Im(RK), (3.15) gives 


V T>u{x)dx 


\K\Xku = \K\XtMu)\k, 


which shows that cj^(P,T, $) = 0, and thus that (2.36) holds. The coercivity, limit-conformity and compactness of 
(Pm)meN therefore follow from Proposition 2.31. Since HMM gradient discretisations are LLE gradient discretisations, 
the consistency of (Pm)mgN readily follows from Proposition 2.14, after noticing that the regularity assumption on 
(Pm)meN gives a bound on (regLLE(®m))meN- □ 
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3.6.2 The SUSHI scheme 

With the notations of Section 2, the SUSHI scheme [33] is the gradient scheme obtained by a barycentric condensation 
of the HMM gradient discretisation. For simplicity, we only consider here the case when all face unknowns are elimi¬ 
nated (the “SUCCES” version in [33]), although more accurate methods could be used in the case of coarse meshes in 
heterogeneous domains. 

1. Let T be a polytopal toolbox of Q in the sense of Dehnition 2.24. We first define an HMM gradient discretisation 
D = (Xx>_o, V-D, Hri), as in the section above, for which I — A4 U £. 

2. We introduce = A4 U fext, and for all a £ fint we select C and introduce barycentric coefficients 
such that 

^ /3f = 1 and PiXi, 

which corresponds to (2.17). 

3. The SUSHI gradient discretisation is the corresponding barycentric condensation of T> in the sense of Dehnition 
2.15. 

4. We have H^^ba = TIti since this reconstruction is only built from the values at the centres of the cells. 

The property (V) for this barycentric condensation of the HMM method is a consequence of the property (V) for the 
HMM method and of Theorem 2.18, assuming that (reg]3,^(Ilm) -I- (^T>m + 9Tm)m£'N is bounded. 


3.7 Nodal mimetic finite difference methods 

The nodal MFD method (nMFD) is described in [14]. We present here a gradient discretisation which enables us to 
write the nMFD method for (2.1) as the gradient scheme (2.3). Let T be a polytopal mesh of Q in the sense of Dehnition 
2.21. For each K £ A4 we choose non-negative weights ((jJxjveVK such that the quadrature 


/ w(x)dx « 2^ ‘^kU’(v) 
Jk 


(3.16) 


veVjc 


is exact for constant functions w, which means that X^vgVj^ = 1 ^1- For each face a £ £k H Tint, we also choose 
non-negative weights (^^(vgv,, such that the quadrature 


/ ui(a:)ds(a;) « 

vGV„ 


(3.17) 


is exact for affine functions w. This is equivalent to X^vev “ I'^l = |cr|®CT. We will assume the 

following property on these weights. This property is not required in the construction of the nMFD, but it is used to 
identify the nMFD with a gradient scheme. We note that this assumption is not very restrictive, since it holds for any 
natural choice of weights for (3.17). 


yK £ AI, Vv £ Vk , da £ £k,v such that / 0, 


(3.18) 


where £k,v = {ct £ £k '■ v £ Vo-} is the set of faces of K that have v as one of their vertices. For each cell K G A4, we 
dehne its centre as 

1 V 

' ' vGVk 


xk 


(3.19) 


and we select a partition (I/ff,v)veVjc such that 


Vv£ Vif, |Vk,v| = ^ 


\D. 


= rf I] 


(3.20) 




We can then again use Definition 2.11 to construct the gradient discretisation D — (X-^.o, Vr>, Hxj) corresponding to the 
nMFD scheme. 

1. The set of geometrical entities attached to the dof is 7 = V and the set of approximation points is S = 7. We set 
7n = V n H and Ign = V n d£l. The partition is 77 = (VK,v)ifeAi,vgVic, and for U = Vk,m we let lu = Vk- 

2. For U = Vk,v we choose in the reconstruction (2.13) the functions 


This leads to 


Va: £ 7/, Vv' £ Vif , Q;v'(a:) := 

Vu £ Xt>,o, V77 £ M, Va; £ K, YVtiv{x) = vk ■= ^ i-^KVv 

' ' vevjf 


(3.21) 

(3.22) 
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3. In a similar way as for the HMM method, the reconstructed gradient is the sum of a constant gradient in each cell 
and of stabilisation terms in each VK,y We set 


VweXc.o, 'iK£M, Vkv 



(3.23) 


and 

Vw € Xt>,o , VVk,v £U , Va: G Vk,v , Qvk ■,v{x) — \/ kv + -—(Qx (n))]vNx,v (3.24) 

’ flK 

where 

• Nx.v = ‘^o-nx.CT, 

• Qk(v) = (Uv - Ux)veVjf with Ux = I^veVx “ (3-22), 

• Rk ■ is the linear mapping described by Rk{0 — (^^,v(?))v6V/f with 

Rx,v(C) = Cv — VxC • (v — xk), (3.25) 

where VxC is defined as in (3.23), and xk is the centre of K defined by (3.19), 

• Ck is an isomorphism of the space Im(Rx) C R^^. 

4. The proof that |lVr> • ||_LP{n)£i is a norm on X-d,o is the consequence of the following inequality, which is obtained 
in the same way as in [29, Lemma 5.3]: defining uk = EveVjf ‘^x^iv (see (3.22)), we have 




Uv — Uk 
hx 




(3.26) 


with Cs only depending on an upper bound on 9j- and on the regularity (jy. This factor is defined as the smallest 
number such that, for all A' G Af and all ^ G R^^, 


ct> 

vSVk 


Rx,v(C) 

Kk 


< E 

vGVjf 


[TxRx(C)]v 

hx 


< Co E 

vSVjf 


Rk,w{0 ’’ 
hx 


(3.27) 


Under Assumption (3.18) it is proved in [27] that the gradient scheme (2.3) obtained from this gradient discretisation is 
identical to the nMFD method of [14] for (2.1). 


Remark 3.9: The second equality in (3.20) comes from \Dx,a\ = and this choice of [Vx,v[ is compatible with the 

requirement that X^veVj^ = 1^1- The detailed construction and geometric properties of Vx,v are not needed for 

the analysis of the method or for its implementation. The only required information is the measure of this set. 

Other choices of Vx.v are possible. For example, we could take all (Vx,v)vgVjc of same measure , and Property 

iV) would still be valid. However, a stronger assumption than (3.18) would be required to ensure the coercivity of the 
corresponding gradient discretisations; we would need with c > 0 not depending on K or v. 


Remark 3.10: Contrary to the HMM gradient discretisation, the nMFD gradient discretisation does not have a piecewise 
constant reconstruction for the natural choice of unknowns, nor for any obvious choice of unknowns. It should therefore 
be modified, e.g. by mass-lumping as in Section 2.4, to be applicable in practice to certain non-linear models. 


Remark 3.11: In the case of octahedral meshes, if the stabilisation term in Vr> is set to 0, then the discrete space and 
gradient of the nMFD gradient discretisation are identical to the discrete space and gradient of the octahedral gradient 
discretisation (and thus of the CeVeFE-DDFV method on degenerate octahedra). The only difference remains in the 
definition of the reconstruction \It>- 


Remark 3.12: In [11], the authors construct a vertex-based compatible discrete operator scheme; they show that it belongs 
to the nMFD family of schemes [11, Section 3.5] and that it satisfies the coercivity, consistency and limit-conformity 
properties [11, Section 4.5]. 


The regularity of sequences of nMFD gradient discretisations (iDm)mGN is defined as the regularity of the underlying 
polytopal meshes (7X)meN (see Definition 2.30) and the boundedness of (Ci5m)meN. 

Proof [Proof of the property {V) for the nMFD gradient discretisation[ As in previous proofs, we drop indices m from 
time to time. We define a control $ of P by T, in the sense of Definition 2.29, by 


1 

W\ 


vGVjf 


and 


Vcr G T, ^(w)ct = 


1 


E 

vev„ 


'iK G M, ^(m)x = ux 


(3.28) 
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Let us prove (2.34). Since cu)) = |o-| we have ^(m)^ — ^{u)k = Ev6v„ ‘^ct(wv — uk)- Therefore, using Jensen’s 

inequality and the fact that ^ we hnd 


MdK,c 

o-GSk 


^{u)a - ^iu)K 


dh 


< ^ dit.a ^ 

(t^Sk v^Vo- 


Uy — UK 


<0^rYl E 

(tG€k vev<7 


Uy — UK 


— E ( E 

veVjf \o-e£jf „ 


Mv — UK 


hf 


= e^rd Y. 1^^.''! 

veVjf 


ttv — Uk 


hi, 


We conclude the proof of (2.34) thanks to (3.26). Since II-dm = n7-4>(u), we have lo^{'D,T, $) = 0 and (2.35) follows. 
For K G M we have Vkm = (Vr^(u))|if- Therefore 


[ VT,u{x)dx^\K\{VT^{u))lK + 
J K 


'Y, J?x(QK(n))lv 'Yj 

crG€x,y 


(3.29) 


Similarly as for the HMM method, for any rj G Im(i?ic) we have X^veVjf X^o-ef^ uJa'n.K,iT = 0. Hence, the last term 
in (3.29) vanishes and (2.36) holds since {'D,T, ^) = 0. Hence the hypotheses of Proposition 2.31 are verified, which 
shows that (Vm)meN is coercive, limit-conforming and compact. 

By noticing that regLLj;(iDm) remains bounded by regularity assumption on (Vm)meN, the consistency of (iDm)meN is 
an immediate consequence of Proposition 2.14 since nMFD gradient discretisations are LLE gradient discretisations. □ 


3.8 Vertex approximate gradient (VAG) methods 


Successive versions of the VAG schemes have been described in several papers [36, 38]. VAG methods stem from the 
idea that it is often computationally efficient to have all unknowns located at the vertices of the mesh, especially with 
tetrahedral meshes (which have much less vertices than cells). It is however known that schemes with degrees of freedom 
at the vertices may lead to unacceptable results for the transport of a species in a heterogeneous domain, in particular 
for coarse meshes (one layer of mesh for one homogeneous layer, for example). The VAG schemes are an answer to this 
conundrum. After all possible local eliminations, the VAG schemes only has vertex unknowns, and it has been shown to 
cure the numerical issues for coarse meshes and heterogeneous media [38, 37, 39]; this is due to a specific mass-lumping 
that spreads the reconstructed function between the centre of the control volumes and the vertices. Let us remark that 
the original version of the VAG scheme in [36] uses the same nodal formalism as Section 3.7, but has been shown in the 
FVCA6 3D Benchmark [42] to be less precise than the version presented here. 

The VAG scheme is defined as a barycentric condensation and mass-lumping of the Pi gradient discretisation on a 
sub-tetrahedral mesh. 


1. Let T = (M,£,P,V) be a polytopal mesh of in the sense of Definition 2.21, except the hypothesis that the 
faces a G £ are planar which is not necessary here. We define a tetrahedral mesh by the following procedure. For 
any K G A4, a G £k, and v, v' G Vo- such that [v, v'] is an edge of a, we define the tetrahedron Tx.o.v.v' by its four 
vertices xk,X cr,y,\i' (see Figure 3.8), where the point Xa corresponding to the face a is 


X(y 


1 

Card(Vo) 


E- 

vGVo 


(3.30) 


We denote by the simplicial mesh corresponding to these Tic.o.v.v'- 

2. We let U = {X^ be the Pi gradient discretisation defined from 7^ as in Section 3.1.1 for k — 1. For 

the gradient discretisation T>, the set I of geometrical entities attached to the dof is / = MUV U £, and the set of 
S of approximation points of is S' = {{xK)KeM, (v)vev, (®<t)o-g£). We define £> as the barycentric condensation 
of £> (see Dehnition 2.15) such that 7®* = A1 U V and the degrees of freedom attached to £ are eliminated by 
setting Her = Vo- and the coefficients — l/Card(VCT) for all v G Vo, which are precisely the coefficients in (3.30). 

3. The VAG scheme is the gradient discretisation V obtained from the gradient discretisation 72^* by performing a 

mass-lumping in the sense of Dehnition 2.19. We split each tetrahedron into three parts crw'y 

and (whose detailed geometry is not needed), and we let Vk be the union of all {T^^r v v')o-,v.v' and K be 

the union of all (Tk o vv')if,o-,v'• This leads to 


Vn G Xx>,o ■ UtiV = E vkXVk + E^''A'4- 

K£M vSV 


The regularity of a sequence of VAG gradient discretisations is dehned as the regularity of the underlying tetrahedral 
meshes in the sense of Dehnition 2.30. We can check that reg0,^(7J ) remains bounded by a non-decreasing function 
of Oj-T, and the proof of the property iV) for VAG gradient discretisation is thus a direct consequence of the results 
in Section 2, especially Theorem 2.18 (properties of the barycentric condensation) and Theorem 2.20 (properties of 
mass-lumped gradient discretisations). 
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Fig. 5: Definition of tetrahedron TK,cr,yy in a mesh cell K. 


4 Conclusion 

We gave here a brief presentation of gradient schemes, a generic framework for the convergence analysis of several 
numerical methods for various diffusion models. This framework is based on the notion of gradient discretisations (a 
triplet of discrete space, reconstructed gradient and reconstructed function) and on core properties they must satisfy to 
ensure the convergence of the corresponding gradient schemes. We provided generic tools to prove that given numerical 
schemes can be associated with gradient discretisations that satisfy the core properties: local linearly exact gradients, 
barycentric condensation, mass lumping and control by polytopal toolboxes. We then showed that several classical 
methods are gradient schemes: conforming Pfc hnite elements (and mass-lumped Pi finite elements), non-conforming Pi 
finite elements (with or without mass-lumping), RT*, mixed finite elements, multi-point flux approximation 0-scheme, 
discrete duality finite volumes, hybrid mimetic mixed methods (which contains hybrid mimetic finite differences), nodal 
mimetic finite differences, and the vertex approximate gradient scheme. All these schemes have been shown to be 
associated with gradient discretisations that satisfy the required core properties. 

Ongoing works concern the adaptation of the gradient scheme framework to some more general operators. In the case 
of the incompressible Stokes equations [26], it has been possible to obtain convergence results which simultaneously hold 
for the Taylor-Hood scheme, the Crouzeix-Raviart scheme and the MAC scheme. Results in this direction have also been 
obtained on the elasticity problems [30]. Some interesting questions remain open even in the case of the Laplace equation. 
For example, it is still not known whether discontinuous Galerkin schemes fall into the gradient scheme framework, i.e. 
if it is possible to construct a gradient that gathers the consistent part of the discontinuous Galerkin gradient and the 
jumps penalisation. As shown in the study of HMM methods[29, Eq. (5.11)], such a construction would require some 
form of orthogonality property between these two components of discontinuous Galerkin schemes; further investigation 
is necessary. 

Acknowledgement Acknowledgement: the authors would like to thank Boris Andreianov for insightful discussions 
during the course of this research. 
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